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^ . ABSTRACT 

Aims. We present the first direct comparison of the distribution of the gas, as traced by the [O i] 6300A emission, and the dust, as 
traced by the 10 /im emission, in the planet-forming region of proto-planetary disks around three intermediate-mass stars: HD 101412, 
HD 135344 B and HD 179218. 

Methods. A'-band visibilities were obtained with VLTI/MIDI. Simple geometrical models are used to compare the dust emission to 
high-resolution optical spectra in the 6300 A [Oi] line of the same targets. 

Results. HD 101412 and HD 135344 B show compact (< 2 AU) 10 /jm emission while the [Oi] brightness profile shows a double 
peaked structure. The inner peak is strongest and is consistent with the location of the dust, the outer peak is fainter and is located at 
^ ' 5-10 AU. In both systems, spatially extended PAH emission is found. HD 179218 shows a double ring-like 10 /jm emission with the 

^ first ring peaking at ~ 1 AU and the second at ~ 20 AU. The [O i] emitting region is more compact, peaking between 3-6 AU. 

. Conclusions. The disks around HD 101412 and HD 135344 B appear strongly flared in the gas, but self-shadowed in the dust beyond 

0^ ' ~ 2 AU. The difference in the gas and dust vertical structure beyond 2 AU might be the first observational evidence of gas-dust 

CO ' decoupling in protoplanetary disks. The disk around HD 179218 is flared in the dust. The 10 /jm emission emerges from the inner rim 

and from the flared surface of the disk at larger radii. No dust emission is detected between ~ 3 - 15 AU. The oxygen emission seems 
also to come from a flared structure, however, the bulk of this emission is produced between ~ 1 - 10 AU. This could indicate a lack 
of gas in the outer disk or could be due to chemical effects which reduce the abundance of OH - the parent molecule of the observed 
[O i] emission - further away from the star. It may also be a contrast effect if the [O i] emission is much stronger in the inner disk. 
We suggest that the three systems, HD 179218, HD 135344 B and HD 101412, may form an evolutionary sequence: the disk initially 
flared becomes flat under the combined action of gas-dust decoupling, grain growth and dust settling. 
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1 . Introduction gas that forms the disk undergoes changes in its composition and 

size. Infrared surveys show that after a mean age of 3 Myr the 

Many pre-main-sequence stars are characterized by excess in- -^^^^^j. disk is cleared of dust. Viscous accretion, photo- 

frared emission above the stellar photospheric level which, de- evaporation and planet formation are the likely mecha nisms re- 

pending on the evolutionary state of the system, may start in the sponsible for this phenomenon (see review in Henning|2008ll. 
near infrared (1-5 fim) or at longer wavelengths. The dust dis- 
tributed around the young star is responsible for this emission. 

The dust particles absorb a large fraction of the short wavelength a circumstellar disk is believed to be the locus where planet 
stellar photons and re-emit them at infrared wavelengths. This formation takes place. The large number of recently discovered 
dust is believed to be confined to a disk-like structure which exo-planetary systems and the variety of such systems raised 
forms in the early phase of star formation as a result of the con- many new questions about the structure and evolution of pre- 
servation of the angular momentum of the parental cloud. Disk toplanetary disks. Of particular interest for the understanding 
formation is followed by a longer phase of disk accretion during of disk evolution and planet formation is the coupling of gas 
which disk material accretes onto the young star at a typical ac- and dust in disks. A main assumption underlying essentially 
cretion rate of 10 ^ - 10 M© yr The interstellar dust and all proto-planetary disk models is that dust and gas are ther- 

mally coupled. It is an open question whether this assumption 

Send offprint requests to: D. Fedele holds on the disk surface and efforts have been made to improve 

* Based on observations collected at the European Southern disk models by takin g into a ccount the dust-gas decoupling (e.g. 

Observatory, Paranal, Chile (Proposal ID: 077.C-052I, 077.C-0263) Kamp & Dullemond |2004| l. While gas and dust are thermally 
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Table 1. Properties of the programme stars. The column "Meeus group" gives the classification m flared (group I) and self-shadowed 
(group II) disks by Meeus et al. (20011. Column "[Oi] extent" reports the extent (minimum and maximum radius) of the [Oi] 
emitting region as derived in paper I. The disk position angle (PA) and inclination are taken from paper I. 



Star 


RA 


DEC 


Sp.T 


M,„,-„. 




Distance 


Meeus Group 


PA 


Inclination 


[Oi] extent 




(J2000) 


(J2000) 




[Mq] 


[Jy] 


[pc] 




n 


n 


[AU] 


HD 101412 


11:39:44.46 


-60:10:27.7 


AOIIIe 


2.3 


3.22 


160 


II 




30 


0.15 - 10 


HD 135344 B 


15:15:48.44 


-37:09:16.0 


F4Ve 


1.7 


1.59 


140 


I 


100 


45 


0.1 - 100 


HD 179218 


19:11:11.25 


+ 15:47:15.6 


B9e 


2.7 


23.4 


240 


I 


10 


40 


0.4 - 50 



coupled in the disk interior, in the low density environments of 
the disk surface layer, the two components may decouple. 

The detailed structure of the disk surface temperature in 
the presence of gas-dust decoupling was studied by Jonkheid 
et al. ((2004), Kamp & Dullemond (|2004] l and Nomura & 
Millar (|2005). Different heating/cooling processes act at differ- 
ent heights in the disk. In particular, very high in the atmosphere, 
with particle densities as low as n < lO^cm"^ (Ay < 10"^ mag), 
the gas temperature is set by the balance of photoelectric heat- 
ing and fine structure line cooling of neutral oxygen (Jonkheid 
et al. 120041 Kamp & Dullemond |2004| l. In the upper layers of 
protoplanetary disks the gas temperature may exceed the dust 
temperature. At small radii (< 50 AU) the temperature of the gas 
above the disk photosphere may reach ~ 10** K. At larger radii 
(> 50 AU) the gas can become as hot as a few hundred Kelvin 
(Kamp & Dullemond |2004| l . However, it is not obvious that such 
a hot and tenuous disk atmosphere can remain as a stable struc- 
ture of a disk. The difference in gas and dust temperature as well 
as other processes (e.g. disk wind, disk evaporation, dust coagu- 
lation and settling) may lead to a physical decoupling of gas and 
dust in disks. 

In this paper we present the first direct comparison of the 
dust and gas emission of three pre-main-sequence stars: HD 
101412, HD 135344 B and HD 179218. These are intermediate- 
mass (1.7-2.7 Mq) stars belonging to both group I (flaring disk: 
HD 135344 B and HD 179218) and group II (flattened di sk: HD 
101412) according to the classification of Meeus et al. (120011 1. 
In a previous paper (van der Plas et al. ,2008 hereafter paper 
I) we presented high resolution spectroscopy {jj - 77000) of 
the optical [Oi] 6300.304A Une with VLT/UVES of these three 
stars. Here we present A^-band interferometric observations ob- 
tained with VLTI/MIDI aimed at spatially resolving the mid- 
infrared emitting region of the disk. We will also compare the 
thermal coupling between dust and gas in the disks suiTound- 
ings our three targets by comparing the MIDI observations to 
a simple phenomenological model derived from the [O i] data. 
The properties of the programme stars are reported in Table [T] 
The paper is organized as follow: section 2 briefly summarizes 
the results of paper I; section 3 contains information about ob- 
servations and data reduction; in section 4 and 5 we report re- 
spectively an analysis of the interferometric measurements and 
the comparison with the optical data of paper I. The discussion 
and conclusion are given in section 6. 

2. [Oi] emission: paper I 

The [Oi] 6300A emission in young intermediate-mass stars is 
thought to be caused by photodissociation of the OH molecules 
present in the protoplanetary disk by stellar UV photons. A frac- 
tion of the resulting excited oxygen atoms is in the upper state 
('D2) of the 6300A transition. The [Oi] 6300A fine is thus non- 
thermal and strongly sensitive to the UV radiation from the cen- 



tral star (Acke et al. '20051. For this reason, the oxygen line traces 
only the surface layers of the disk that are directly exposed to the 
stellar radiation field. 

The Doppler broadened oxygen emission profile was trans- 
lated into an amount of emission as a function of distance from 
the central star assuming Keplerian rotation of the gas. The re- 
sulting radial profile of the [O i] emission is in agreement with 
the expected disk shapes as derived from their spectral energy 
distribution (SED) according to the phenomenological classifi- 
cation of Meeus et al. (12001I) . For all targets the oxygen emis- 
sion starts at velocities corresponding to their dust sublimation 
radius and extends up to radii of 10 - 90 AU. The [Oi] radial 
profile shows a double peak structure in the case of HD 101412 
and HD 135344 B, with a first stronger peak at ~ 1 AU, and a 
second weaker peak at ~ 5-10 AU. In HD 179218 the [O i] has a 
single peak at ~ 3 - 6 AU. The normalized intensity [O i] radial 
profiles for the three stars are shown in Figs.[8ll9land[T0l 

3. Observations and data reduction 

MIDfl (Leinert et al. |2003] l is the mid-infrared (8-13 yum) 
beam combiner of the VLT interferometer on Cerro Paranal, 
Chile. MIDI can combine the light of any pair of the 4 ESO 
VLT telescopes (UT, 8.2-m). MIDI allows spectrally resolved 
observations, that is, the signal is dispersed by either a prism 
(A/1//1 X! 30) or a grism {AAjA ^ 300). For our project we used 
the low resolution (prism) mode and two nearly perpendicular 
basehnes: UT1-UT2 {B = 57m, PA = 26°) and UT3-UT4 {B = 
62m, PA = 111°). Taking advantage of the earth rotation (which 
modifies the projection of the baseline on the sky) we observed 
our targets at different sidereal times, to optimize the coverage 
of the (u,v) plane. 

Two modes are available for the acquisition of the interfero- 
metric signal: "High Sens" and "Sci Phot". With the first mode, 
the photometric signal (i.e. the pure, not combined, 8-13 jum 
spectrum of the source coming from the two telescopes) is ac- 
quired soon after the interferometric signal. In Sci Phot mode, 
the two signals are acquired at the same time. We used both 
modes depending on the N-band luminosity of the source. The 
log of MIDI observations is reported in Table |2l The data were 
reduced with the standard data reduction software MIA-hEWS 
vl. 5. fusing the Expert Work Station (EWS). Detailed doc- 
umentation on the reduction procedure can be found in Ratzka 
(120051) and Jaffe (120041 1 or in the cited web pages. 

Standard stars (HD 107446, HD 129456 and HD 188512) 
were observed close in time and in airmass to the science tar- 
gets. The angular diameter of these stars is well known and they 
can be used as calibrators for the transfer function of the instru- 
ment. Other calibrators observed during these nights were used 

' www.eso.org/instruments/midi 

- www.mpia-hd.mpg.de/MIDISOFT 

3 www.strw.leidenuniv.nl/ koehler/MIA+EWS-Manual 
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Fig. 1. Transfer function of VLTI/MIDI. For each night, the different lines correspond to different calibrators observed with the 
prism in both modes: HS (solid line) and SP (dashed line). The overall shape of the transfer function, for a given mode, does not 
drastically change during the night. The vertical grey regions are characterized by strong atmospheric absorption bands. Residuals 
of the data reduction may be present within these regions. 



to check the stability of the transfer function. The error on the 
calibrated visibility of the science targets is the standard devia- 
tion resulting from calibration with different calibrator stars for 
each night. Fig.[T]shows the transfer function (i.e. the instrumen- 
tal visibility of the calibrators) for the six nights of observations. 
For each night, all the calibrators observed with the prism in 
both modes (HS and SP) are plotted. The overall shape of the 
transfer function with wavelength varies by less than 10% dur- 
ing the night. The standard stars were also used to flux calibrate 
the mid-infrared spectrum of the science targets. Conversion fac- 
tors from counts to Jy were computed from theoretical spectral 
energy distributions (SED). Energy distributions were derived 
by matching the spectral type of the calibrators to stars in the 
Cohen list (Cohen et al. J999 ) and scaling the spectrum with the 
ratio of their IRAS 12 jum fluxes. 

The averaged MIDI spectra of the three targets are shown 
in Fig.|2l The Spitzer spectrum is plotted for comparison. Each 
spectrum was scaled in order to match the Spitzer spectrum at 
9 jum. MIDI and Spitzer spectra agree well. A systematic devia- 
tion of the MIDI spectrum from the Spitzer one is visible at the 
red edge of the N-band. This is due to a non-perfect correction 
of the MIDI spectrum for atmospheric absorption. 



3A. Note on the data reduction of HD 135344 B 

The data reduction process of HD 135344 B deserves a more 
detailed description. The brightness of this target is close to the 
sensitivity limit of the instrument. To check the quality of the 
data we used both the coherent (EWS) and incoherent (MIA) 
method. Given the faintness of the source, we used a higher 
threshold for the "good scans" in MIA (60 - 70 %). In all cases, 
small differences at the level of 10 - 20 % in the calibrated vis- 
ibility are present. In particular, on the night of 2006 July 12, 
EWS yields a lower visibility at short wavelengths. After dif- 
ferent checks, the problem seems to be the fixed mask used to 
extract the signals with EWS. Both the bi-dimensional interfer- 
ometric and photometric signal are slightly offset compared to 
the position of the EWS mask. For this reason we preferred to 
adopt the mask computed by MIA. With this procedure the two 
data reduction packages produce very similar results which are 
consistent with a precision of < 10%. 

4. Data analysis 

The quantity measured by an interferometer is the Fourier trans- 
form of the brightness distribution, i.e. the complex visibility. 
Due to atmospheric turbulence, however, the phase of this com- 
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Table 2. Summary of VLTI/MIDI observations. Column "Mode" indicates the beam combination mode of the interferometric signal: 
High-Sens (HS) or Sci-Phot (SP). Columns "B^" and "PA" list the length of the projected baseline and its position angle (east of 
north). 



Night 


UT 




Mode 


A/AA Baseline 


fix 


PA 


Airmass 


Calibrator 




[hh:mm] 






[m] 


n 






HD 101412 


2006-05-15 


23:38 


- 00:02 


SP 


30 UT3-UT4 


59 


102 


1.2 


HD 107446 


2006-05-16 


01:55 


-02:18 


SP 


30 UT3-UT4 


62 


110 


1.3 


HD 107446 


2006-06-13 


00:18 


- 00:34 


SP 


30 UT1-UT2 


43 


30 


1.3 


HD 107446 


HD 135344 B 


2006-05-13 


04:32 


-05:15 


HS 


30 UT3-UT4 


62 


112 


1.0 


HD 129456 


2006-06-13 


23:26 


- 23:52 


HS 


30 UT1-UT2 


56 


27 


1.3 


HD 129456 


2006-07-12 


03:14 


- 03:47 


HS 


30 UT1-UT2 


47 


10 


1.3 


HD 129456 


HD 179218 


2006-05-16 


05:25 


- 05:44 


SP 


30 UT3-UT4 


58 


128 


1.7 


HD 188512 


2006-06-13 


06:42 


- 06:56 


SP 


30 UT1-UT2 


48 


42 


1.3 


HD 188512 


2006-06-13 


08:07 


- 08:25 


SP 


30 UT1-UT2 


53 


40 


1.5 


HD 188512 


2006-07-08 


06:45 


- 07:03 


SP 


30 UT3-UT4 


44 


84 


1.6 


HD 188512 



plex number is corrupted and cannot be retrieved. The amplitude 
of the complex visibility - the visibility - can be measured and is 
compared to the corresponding quantity of theoretical brightness 
distributions. The visibility is a two-dimensional function in the 
so-called u, v-plane, where u and v are the spatial frequencies 
B/A(B is the baseline length between the two telescopes, and A 
the wavelength) in two perpendicular directions. In the case of 
VLTI/MIDI, a single observation covers a range of spatial fre- 
quencies due to the spectral capability of the beam combiner. 
However, one should keep in mind that the brightness distribu- 
tion of the target may vary from wavelength to wavelength. 

Figs. |3] m & |6] show the calibrated MIDI visibilities of HD 
101412, HD 135344 B and HD 179218 respectively. We plot 
the visibility versus spatial frequency rather than wavelength, 
in order to make a direct comparison with the spatial distribu- 
tion of the [O i] emission measured at 6300A. For convenience, 
we plot the MIDI wavelength scale at the top of the figure. The 
slope of the visibility curve in the A^-band depends on: 1) the 
increase in spatial resolution with decreasing wavelength and 2) 
the wavelength-dependence of the brightness distribution. 

The [O i] visibility curves (dashed line) are calculated from 
the radial intensity profile presented in paper I and are discussed 
later in the text. The comparison of MIDI observations and [O i] 
visibility is presented in Sec. |5] In this section we use a simple 
geometrical model aimed at deriving the geometrical properties 
(inner and outer radius, inclination and position angle) of the 
(mid-infrared) dust emitting region in the three protoplanetary 
disks. 

4. 1 . Uniform ring model 

In a first approximation the brightness distribution of a proto- 
planetary disk might be described by a uniform bright ring (UR). 
The formalism of the UR model is described in Appendix B. We 
emphasize that the UR model is a simple geometrical approxi- 
mation of the dust emitting region. This model does not take into 
account the physical properties of the circumstellar disk (e.g. 
dust sublimation radius, disk outer radius, temperature gradient) 
or the dust properties (e.g. grain size). Such a model has been 
extensively used as an analytic model for near-infrared interfer- 
ometric observations of disks around young stars (e.g. Millan- 
Gabet et al. 2001 Eisner et al. l2004l |20^ as well as mid- 
infrared interferometric data of planetary nebulae (Chesneau et 



al. |2006l ). The MIDI visibilities were fitted using an inclined uni- 
form ring (UR) model. 

We have searched for the best set of model parameters by 
minimizing the of the visibility measurements. The values 
of the best fit parameters and associated errors were computed 
using a Monte Carlo simulation. Assuming a normal error dis- 
tribution of the measured visibility, we simulated 100 random 
data-sets around the observed values. The best fit parameters and 
errors correspond to the mean and the standard deviation of the 
100 fits. 



4.2. HD 101412 

HD 101412 (Fig. [3]) is barely resolved at all epochs. Along 
the three baseline position angles the visibility is lower at 
short wavelengths and increases from 8 fim to 9 fim - 10 //m. 
Diff'erences in the absolute value and shape in the three mea- 
surements may suggest an asymmetric emission. In particular, 
HD 101412 is more resolved (i.e. shows a more extended emis- 
sion) with the shortest baseline, 43 m at 30°. Along this direction 
the visibility goes from 0.55 at 8 fim to 0.9 at 13 fim. Around 
11.3 jim the disk appears larger than at adjacent wavelengths 
(with a small dip in visibility). This coincides with PAH emis- 
sion. The high visibility values suggest a compact (mid-infrared) 
dust emission. The decreasing visibility with decreasing wave- 
length is due to the increase in spatial resolution. The three visi- 
bility measurements of HD 101412 are best fitted by a UR model 
(Fig. [3] solid line) of the inner and outer radius of 0.4 and 1.9 
AU, an inclination of 80° and position angle of the disk major 
axis of 38°. The almost edge-on orientation found here explains 
why the [Oi] line profile shows structures that are not seen in 
other [Oi] profiles of Herbig stars; thanks to the high inclina- 
tion, the projected velocities are close to the real velocities, and 
the Doppler-induced spread in the spectrum is optimal. 

According to the classification of Meeus et al. (I2001I I. the 
SED of HD 101412 suggests a self-shadowed disk geometry. In 
this case the mid-infrared disk emission is expected to arise only 
from the inner rim of the disk. This is consistent with the small 
size measured with MIDI. 
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Fig. 2. MIDI spectra (diamonds) of HD 101412, HD 135344 B 
and HD 179218. The MIDI spectra were scaled to match the 
Spitzer IRS spectra (red Une) at 9 yum. 



4.3. HD 135344B 

HD 135344 B is resolved with all baselines (Fig.gli. Along the 
47 m baseline (position angle 10°) the visibility decreases with 
wavelength. Along the second and third baselines, 56 m and 
62 m respectively, the visibility shows a similar shape: nearly 
constant between 8 //m - 11 yum with a small bump around 
9.2 - 9.4 yum, a drop at 11.3 //m and almost constant at longer 
wavelength. Despite the higher resolution, the source is less re- 
solved at short wavelengths. This means that at longer wave- 
lengths the target is significantly more extended due to emission 
of the colder outer disk which only starts to show up at these 
wavelength. This specific pattern is sometimes observed in mid- 
infrared interferometric observations of young PMS stars (e.g. 
VI 647 Orionis, Abraham et al. i2006i FU Orionis, Quanz et al. 
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Fig. 3. Calibrated visibility amplitude at different baselines for 
HD 101412. The dashed line in the visibility panel represents 
the visibility computed from the intensity versus radius profile of 
the [O i]6300A emission line. The solid line is the best fit of the 
inclined UR model. Visibility errors are computed as standard 
deviations resulting from calibrations with different calibrator 
stars. 
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Fig. 4. Same as Fig.|3]for HD 135344 B. 



I2OO6I 1. Reasonable fitting of this visibility shape is obtained with 
a (physical) Keplerian, flared disk model where temperature and 
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Fig. 5. Correlated spectrum of HD 135344 B measured on 2006 
May 13 (dotted), 2006 June 13 (dashed) and 2006 July 12 (dot- 
dashed). The spectra were shifted for clarity. The mean, single 
telescope, spectrum is also plotted (solid line). The absence of 
the PAH emission feature at 1 1.3 fim in the correlated spectrum 
is a signature of spatially extended PAH emission compared to 
the continuum emission. 

surface density are prescribed by a broken power law distribution 
(T oc r^i\ 2 oc r^P), i.e. the disk's surface is cooler at larger radii. 
The result of the UR model (Fig.lH solid line and Table O are: 
inner radius 0.05 AU, outer radius 1.8 AU (assuming a distance 
of 140 pc van Boekel et al. 12005), inclination 60 ° and position 
angle of the disk major axis 180°. The inclination found here is 
slightly higher than the result of Doucet et al. (2 006) (4 5°), ob- 
tained with direct imaging at 20.5 jjm. Dent et al. ( 120051 1. model- 
ing the J=3-2 transition of '^CO, find a smaller disk inclination 
(11°, disk nearly face-on). Nevertheless, the UR model is not 
able to reproduce the decreasing visibility observed with MIDI 
and the disk parameters are not well constrained. 

We converted the MIDI visibility to the Gaussian full width 
half maximum (FWHM) at each spatial frequency. The FWHM 
represents the extent of the emitting region and is defined as: 



IniVisf)) 
' -3.56i/2 



(1) 



where V(sf) is the MIDI visibility at spatial frequency sf. 
A Gaussian distribution is only a crude approximation of the 
real brightness distribution valid at high visibilities. Similarly 
to Quanz et al. (2006), we also fitted the FWHM measured on 
each baseline with an ellipse wavelength by wavelength. The 
best-fit ellipse parameters at three reference wavelengths (9 //m, 
11 fim, 12.5 fim) are listed in Table |4| As expected, the disk 
looks larger at longer wavelengths. The ellipse semi-major axis 
increases from 1.1 AU at 9 fj.m to 1.9 AU at 12.5 jum. The incli- 
nation and position angle varies by ~ 5° and ~ 9°respectively 
from 9 yum to 1 1 jum. The 1 1 fj.m and 12.5 /vm ellipses have simi- 
lar (within ~ 3°) inclinations and position angles. The geometry 
of the mid-infrared dust emission derived with the UR model 
and the ellipse fitting are consistent with each other 

The SED of HD 135344 B is characterized by a relatively 
low infrared excess over the stellar photosphere up to 13 //m 
and by a strong infrared excess beyond 14 //m. Moreover, the 
mid-infrared spectrum (Fig. |2]) lacks silicate features and only 
weak PAH emission is detected. Such an SED is peculiar among 
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Fig. 6. Same as Fig.|3]for HD 179218. The continuous line here represents the best fit two-component disk model. 



group I sources and suggests a more complex disk structure. 
The missing near- to mid-infrared excess and the fast increase 
at longer wavelengths have been attributed to the presence of a 
large dust-free gap by Brown et al. ( 120071 ). They proposed a so- 
called "cold disk" structure for the material around HD 135344 
B. In their model the disk is devoid of dust between 0.45 AU 
- 45 AU (in radius). According to Brown et al., warm dust is 
present in an inner ring between 0. 18 AU - 0.45 AU. The MIDI 
observations presented here are consistent with this model. The 
Gaussian FWHM, if scaled to the distance adopted by Brown et 
al., yields an emitting region between 0.4 - 0.7 AU. Although 
these radii are larger than those found by Brown et al. (note 
that we used a simple Gaussian model), the two estimates agree 
within the uncertainty of the data and the analysis. 

Because MIDI is only sensitive to dust emitting at 10 fim, 
our data set does not allow us to claim the presence or absence of 
cold matter beyond the inner warm disk. The shadow of the inner 
disk may obscure a large portion of the disk which re-emerges 
only at large radii. Both the SED and MIDI observations are only 



sensitive to the dust emitting in the mid-infrared, both the SED 
and our MIDI observations are not able to distinguish between a 
(partially) self-shadowed disk and a disk with a gap. 

The correlated spectrum of HD 135344 B for the different 
measurements is plotted in Fig. |5] At all epochs the correlated 
spectrum mimics the total (single-telescope) spectrum, although 
deviations are visible. The 11.3 //m PAH emission detected in 
the total spectrum is not present in any of the correlated spectra 
which indicates that the PAH 11.3 /vm emission region is com- 
pletely resolved at the angular resolution of our measurements. 

4.4. HD 179218 

Despite being the most distant object of our sample, HD 179218 
is the most spatially resolved source of the three objects. The 
visibility varies with the position angle of the baseline and is al- 
ways below 0.4. Along the first baselines the visibility increases 
rapidly from 8 to 9 fim and decreases afterwards. Along the sec- 
ond and third baselines (48 and 53 m), the visibility is much 
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lower and shows a characteristic sinusoidal modulation. In the 
last case, the visibility increases from a very low value, 0.05 at 
8 fim to 0.27 at 11. 7 fim and decreases afterwards. 

A sinusoidal pattern of visibility is produced if the brightness 
distribution of the target displays sharp edges, e.g. in the case 
of a ring, uniform disk or a binary with two unresolved compo- 
nents. A smooth emission region, such as a Gaussian disk, would 
not produce such variations. We can immediately exclude the bi- 
nary hypothesis: the MIDI observations of HD 179218 cover a 
range in position angles over 90° and the visibility always shows 
a sinusoidal pattern. In the case of a binary, the pattern is PA de- 
pendent and completely disappears when the PA of the baseline 
and that of the binary axis are aligned. We attempted to fit the 
MIDI visibility with the UR model but although this model is 
able to reproduce the sinusoidal modulation of the visibility, it 
fails to reproduce the visibility at the position of the local "min- 
ima". The minimum in the ring model is zero while the observed 
visibility reaches a minimum value of V = 0.05. If these were 
"true" zeroes of the visibility, i.e. the location between two ad- 
jacent lobes of the visibility curve, we would have measured a 
contrast of V ~ 0.007(0. We then attempted to fit the data with a 
two-component model. The first component represents the disk 
emission at small radii and the second component represents 
the disk surface at larger radii. Our simulations showed that the 
first component is only barely resolved. We used a uniform ring 
for both the internal and external component. Using a different 
brightness distribution for the internal component does not im- 
prove the model fit. The fractional flux contribution of the first 
component, //„,, was left as a free parameter in the fit. For sim- 
plicity, fi„,, was assumed to be constant at all wavelengths. The 
best fit parameters are listed in Table |5] Inclination and posi- 
tion angle are in agreement with those found with the single ring 
model (i - 57°, PA - 198°). The inner component contributes 
20% of the total mid-infrared flux. According to this model, the 
disk around HD 179218 is either truncated or obscured by its 
own shadow between 3.2 and 14.5 AU radii from the star (Table 

Previous infrared interferometric observations of HD 
179218 have been carried out by by Leinert et al. (|2004| l and Liu 
et al. (12007t in the mid-infrared and by Monnier et al. (120061) in 
the H-band. Leinert et al. also measured low visibilities with an 
80 meter baseline. Monnier et al. ( |20061 l measured a squared H- 
band visibility of the order of 40%, corresponding to a Gaussian 
FWHM of 3.8 AU. This is in good agreement with our two- 
component model, where the inner component extends up to 3.5 
AU from the star. Liu et al. modeled MMT nulling interferomet- 
ric data by means of a ring of diameter 27 + 5 AU. This large 
ring is, within the uncertainties, consistent with the outer ring of 
our two-component model. 

The agreement between the MIDI data and the simple geo- 
metrical model proposed here is quite good. However, deviations 
from the model are clearly visible. Similarly to HD 135344 B, 
our observations are not able to distinguish between a real gap 
(i.e. a region free of material) and a shadowed region. 

4.4.1. Differential pliase 

As a byproduct of the EWS data reduction software, the relative 
phase between the different spectral channels of MIDI can be 
retrieved. This so-called differential phase is the real differential 

This value, given the finite band-width of the MIDI wavelength 
bins, corresponds to the visibility null reachable by the interferometer 
at low spectral resolution (PRISM). 



phase of the target on the sky, modulated by the atmosphere. The 
largest effect of the latter is an offset of the differential phases, 
which is linear with respect to wavenumber (oc A^^). The EWS 
software subtracts the best-fit linear component of the measured 
differential phase, as this information is lost anyway. Secondary 
atmospheric effects remain: the two light beams, coming from 
the two telescopes, pass through different amounts of air and wa- 
ter vapor before they reach the MIDI beam combiner. However, 
these smooth variations in the measured differential phase are of 
the order of a few degrees. For all targets the differential phase 
is nearly constant with wavelength and centered at zero degrees. 
The only exception is the measurements of HD 179218 on 2006 
May 16 (Fig.|7]). Along this baseline, a large phase variation is 
detected at short wavelengths. The phase varies by 55° between 
8 and 9 fim. A flip of 180° is expected when the visibility ampli- 
tude crosses a null, i.e. when the spatial frequencies covered by 
the observations are located between two adjacent lobes of the 
visibility curve. In the case of HD 179218 the phase variation 
is smooth with lower amplitude. The phase variation is then as- 
cribed to a shift of the photo-center of the emitting source. This 
suggests a deviation from centro-symmetric emission along this 
direction. 

The differential phase plotted in Fig.|2]was not calibrated for 
the effect of water vapor that is present in the delay line tunnel 
of the VLT interferometer The two light beams most likely pass 
through a different amount of air and water vapor before reach- 
ing the MIDI beam combiner. This difference introduces an in- 
strumental differential phase. Fig. [T] also shows the differential 
phase measured with four different calibrators observed during 
the same night with the same acquisition mode (SCI PHOT). 
Since calibrators are (usually) point-like their differential phase 
should be zero. The deviation from zero degree in this case is 
produced by the water vapor The effect of water vapor is present, 
but is very small compared to the phase variation observed to- 
ward HD 179218. 

The variation in the differential phase along the 58 m base- 
line is caused by asymmetric emission. This is in agreement 
with the inclination of the disk (57°) and its flared structure. 
Due to the flaring, the far side of the disk emits more light in 
the direction of the observer. This means that the photo-center 
of this emission (described here by the ring) is no longer co- 
incident with the photo-center of the almost unresolved inner 
disk. It is difficult to derive a quantitative estimate for this off- 
set, mainly due to the loss of the linear differential phase infor- 
mation. Estimating from the magnitude of the differential phase 
variation (~60°), it must be around 4+2 mas. 

5. Comparison to the optical gas emitting region 

In paper I very high resolution optical spectra of HD 101412, HD 
134355B and HD 179218 were presented. In all three targets, we 
spectrally resolve the [O i] 6300A emission line. Assuming that 
the disk is in Keplerian rotation, the [O i] emission line profile 
can be translated into a radial profile, i.e. the intensity-versus- 
radius curve of the oxygen line emission, I[Oi]{R). Figs. [8] |9] and 
[TOl show the normalized [O i] radial profile derived in paper I: 
/[O :](/?) X S(R), where S{R) = 7r(/?i„ - Rfj and and R,, 
outer and inner radius of the ring at radius R. I[Oij{R) x S{R) is 
the total luminosity in a ring at distance R from the central star. 
The [O i] radial profile shown here is based on the average of the 
spectra presented in Paper I. In the case of HD 101412, /[Oi](^) 
was recomputed using a disk inclination of 80° to be consistent 
with the result of Sec. 5. Using a larger inclination has the effect 
of shifting the emission further away from the star as can be seen 



Fedele et al.: The structure of protoplanetary disks 



9 




Fig. 7. DiflFerential phase of HD 179218 measured on 2006 May 
16. A 55° phase variation is detected at short wavelengths. The 
four continuous lines are the differential phases for four different 
calibrators observed during the same night. The plotted phases 
are not corrected for the effect of water vapor in the delay lines 
of the VLTl. See Sec. 4.4.1 for details. 

Table 3. Best fit parameters of the UR model for HD 101412 and 
HD 135344 B. Ri„ and Rout are the inner and outer ring radii and 
are given in milliarcsec (first row) and AU (second row). The 
last column lists the value of the reduced;^'^. Values close to one 
indicate a good agreement between model and observations. 



Table 4. Best-fit ellipse parameters for HD 135344 B. The el- 
lipse is fitted to the FWHM measured on the three baselines. 
Ellipse axis are computed assuming a distance of 140 pc. 



Wavelength 


9 yum 


1 1 yum 


12.5 fim 


Semimajor axis [AU] 


1.1 


1.5 


1.9 


Semiminor axis [AU] 


0.6 


0.8 


0.9 


Inclination [°] 


53 


58 


61 


Position angle [°] 


171 


163 


163 



comparing Fig.[8]with Fig. 1 1 of paper 1. The filled region repre- 
sents the position of the dust emission region as determined by 
MIDI. The striped regions indicate the uncertainty on the latter 

In order to compare the [Oi] result with the dust one we 
need to provide a visibility model for the MIDI observations. 
From I[Oi](R), it is possible to derive a brightness distribution 
(I[0i](a,/3), with a and /3 the angular distance from the central 
star) assuming that the line emission from the disk surface is 
axis symmetric. Using I[0i](a,j3) as brightness distribution, the 
visibility model (V[Oi](m,v)) was computed. The expected vis- 
ibility at the spatial frequencies mapped with MIDI was com- 
puted interpolating V[Oi]iu,v) at the corresponding positions in 
the uv-plane. The dashed line in Figs.|3]|4]&|6]is V[Oi](u,v) at 
the spatial frequency covered by MIDI observations. 



Table 5. Best fit parameters and standard deviation of the two- 
component disk model for HD 179218. Radii are given in AU 
assuming a distance of 240 pc. The inclination and position an- 
gle of the two components were assumed to be equal for both 
structures. //„ is the fractional contribution of the inner compo- 
nent to the total flux in the mid-infrared. 



HD 179218 



1" component: 
1" component: 
2"'' component: 
2'"* component: 
inclination (°) 
PAD 

_A 



(AU) 
'ou, (AU) 
J,„ (AU) 



Rour (AU) 



0.3 ± 
3.2 ± 
14.5 
22.6 



0.1 
0.1 
tO.5 
fcO.4 



57 ± 2 
23 + 3 
0.2 ± 0.02 



Star 


Rill 




i 





X' 




[mas/AU] 


[mas/AU] 


[°] 


[°] 




HD 101412 


2.5 ± 0.6 


11.9 ±0.6 


80 ±7 


38 ±5 


OA 




04 ±0.1 


1.9 ±0.1 








HD 135344 


0.35 + 1.79 


12.8 ± 1.4 


60 ± 10 


180 ± 60 


3,, 




0.05 + 0.25 


1.8 ±0.2 









5.1. HD 101412 

As shown in Fig.|3] V[Oi](m, v) is similar to the observed visibility 
along the first baseline (43 m, 30°) and is clearly lower than 
along the other two directions. This is clearer when looking at 
Fig. |8] The dust emission is located within ~ 2 AU, while the 
gas emission appears more extended with a bright peak at ~ 5 - 
6 AU. In terms of radial extent, this means that the [O i] emission 
extends to larger radii than the mid-infrared emission. Beyond 
the first peak the [O i] intensity drops off as one would expect 
for a self-shadowed disk. 



5.2. HD 135344 B 

Similarly to HD 101412 the [Oi] brightness profile of HD 
135344 B shows a double peak with a first stronger peak at small 
radii (~ 0.5 AU) and a second fainter peak at ~ 5-10 AU (Fig. 
|9]l. On the other hand, the 10 fim dust emission seems to come 
from the inner 1-2 AU of the disk surface. However, the [Oi] 
result suffers from overlapping photospheric absorption and thus 
have to be interpreted with more caution (paper I). The diagonal 
striped region in Fig.|9]represents the cold dust detected at longer 
wavelengths by Doucet et al. (|2006l l and by Brown et al. ( 120071 1. 

5.3. HD 179218 

As shown in Fig. |6] the [O i] model is not able to reproduce 
the observed visibility of HD 179218. The observed visibility 
is in all cases lower than Vio ij . Moreover, the oxygen model in- 
creases with wavelength while the observed visibility shows a 
sinusoidal variation. The [O i] intensity peaks at 4 - 6 AU (Fig. 
[Tol l and decreases outwards. According to our two-component 
disk model, emitting dust is present within the oxygen peak. A 
further contribution to the mid-infrared emission, 80% in terms 
of flux, comes from a ring located between ~ 15-23 AU. At this 
distance, the [O i] line intensity is reduced. Our result is consis- 
tent with a flared geometry for the disk around HD 179218; the 
inner, curved, rim is represented by the inner UR while the outer 
UR may represents the flared disk surface at larger radii. The two 
dust emission regions are separated by a "dark" (non emitting) 
region. 



6. Discussion and conclusion 

In this paper we presented the first direct comparison of gas and 
dust emission from the surface of three protoplanetary disks. The 
comparison of high optical spectroscopy with infrared interfero- 
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R[AU] 

Fig. 8. Normalized [O i] 6300A intensity profile for HD 101412. 
Aoi](^) X S{R) is the total [Oi] line luminosity in a ring at dis- 
tance R from the star. This was obtained averaging the results 
of observation Bl and B2 of paper 1. In the specific case of HD 
101412, /[Oi](^) was recomputed with an inclination of 80°. The 
filled region represents the position of the dust emission region 
as determined by MIDI. The striped regions indicate the uncer- 
tainty on the latter. 
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Fig. 9. Same as Fig.[8]for HD 135344 B. The diagonal striped 
region represents the cold dust in a flared geometry detected by 
Doucet et al. (i2006l l and Brown et al. ( 2007 ). This cannot be seen 
with MIDI because it radiates at wavelengths longer than those 
traced with MIDI. 



metric observations gives some insight into the relative size and 
shape of the gas and the dust emitting region. 

For HD 179218 the new MIDI results presented here and the 
UVES results of paper I are in good agreement with a flared disk 
structure. The 10 /im emission comes from two separate regions: 
an inner part located between ~ 0.3 - 3 AU (likely the disk inner 
rim) and an outer part between ~ 15-23 AU (flared disk sur- 
face). No dust emission is detected between 3-15 AU. The [O i] 
brightness profile shows instead a single strong peak centered at 
~ 3 - 6 AU. The difference between gas and dust emitting re- 
gions for HD 179218 may reflects a real different structure for 
the two components or might be caused by other effects such as: 
1) the action of chemical processes that reduce the abundance of 




R[AU] 

Fig. 10. Same as Fig.[8]forHD 179218. 

OH in the outer part of the disk (in this case plenty of gas may 
still exist but we are not able to see it); 2) a contrast effect of 
the [O i] emission. In the latter case the [O i] emission is much 
stronger in the inner disk and it may outshine the [O i] from the 
outer disk. The presence of a dust gap needs further investigation 
as it may indicate possible ongoing planet formation. 

HD 135344 B and HD101412 show a compact 10 //m emit- 
ting region. The oxygen instead has a double peak brightness ra- 
dial profile with the first peak coincident with the dust emission 
and a second peak further away from the star at ~ 5 - 10 AU. 
In the case of HD 135344 B extended dust emission is found at 
20.5 yum (Doucet etal.|2006] l. Modefling the SED of HD 135344 
B, Brown et al. ( I2007I I propose a "cold-disk" geometry where 
the disk is mostly void of dust between 0.4 - 45 AU. Dust is 
present within the inner 0.4 AU and outside 45 AU. This is in 
agreement with our result in the sense that we detect a compact 
mid-infrared dust emission; the dust at larger radii is cold and 
radiates at wavelengths longer than those traced by our observa- 
tions. We suggest that the second [Oi] peak arises either from a 
dust gap or from a dust-free layer above the shadow of the inner 
rim. 

Following the phenomenological classification of Meeus et 
al. C200 1). HD 101412 is a group II source, i.e. the disk is flat 
and self-shadowed in the dust. At the dust sublimation radius 
the disk is puffed-up and shadows the outer disk, which is not 
reached by the stellar UV radiation. The second peak of the [O i] 
emission (Fig. 8) is not expected for such a self-shadowed ge- 
ometry. We suggest that this is the result of a flat disk geometry 
for the dust and a flared geometry for the gas. This might be the 
signature of a different scale height, and hence vertical structure, 
between gas and dust beyond the inner rim. This means that - in 
this object - gas and dust are physically decoupled in the surface 
layers of the disk allowing the gas to emerge from the shadow 
of the inner rim. If this is confirmed, it may have important con- 
sequences for the disk evolution: once gas and dust decouple 
the dust grains settle faster towards the mid-plane. The dust set- 
tling favors the physical separation between gas and dust until all 
the dust settles towards the mid-plane (Dullemond & Dominik 
2004) eventually leaving a region of gas in the upper layers of 
the disk. We thus suggest that HD 179218, HD 135344 B and 
HD 101412 form an evolutionary sequence where the disk, ini- 
tially flared, becomes flat under the combined action of gas-dust 
decoupling, grain growth and dust settling. Interestingly, Acke 
et al. (I2004I I already found that cold grains in the mid-plane of 
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Fig. 11. Simplified sketch of the disk around the three program 
stars. The regions where the dust and [Oi] emissions are ex- 
pected to arise are indicated. The inclination of the disk to the 
line of sight is shown by the dashed arrow at the upper-left cor- 
ner. HD 179218: the disk is large and flared. No dust emission 
is detected between 3-15 AU. Most of the silicate mid-infrared 
emission arises from the warm dust located in the surface lay- 
ers of the outer part while the gas emission is mainly confined 
within 10 AU to the star HD 135344 B: the 10 fim emission 
arises from a compact region close to the star Cold dust is de- 
tected at longer wavelengths. [Oi] 6300A emission is detected 
within the 10 /urn emitting region and further away at 5 - 10 
AU where no dust emission is detected. HD101412: a compact 
dust ring that directly faces the star is responsible for the mid- 
infrared excess. The [Oi] 6300A line arises from the same re- 
gion. Beyond this ring, the disk is shadowed. A dust-free flaring 
layer emerges from the shadow at larger radii. In all three cases, 
gas is likely present within the dust sublimation radius. 



the disk have grown to considerably larger sizes in Meeus et al., 
group II sources than in group I sources, suggesting that disks 
may evolve from a flared to a self-shadowed geometry. Similarly, 
Bouwman et al. (2008 ) find a correlation between the strength of 
the amorphous silicate feature and the shape of the SED, consis- 
tent with the settling of dust as a consequence of grain growth. 
What we add to this picture is a completely independent piece 
of evidence supporting this and perhaps the identification of the 
driving mechanism of this process. 



Appendix A: Extended PAH emission 

Emission from polycyclic aromatic hydrocarbons (PAHs) is typ- 
ically detected towards disks around Herbig stars (e.g. Meeus 
et al. 2001. Acke & van den Ancker 2004; Sloan et al. 2005 & 
Geers et al. |2007) and in some T Tauri stars (e.g. Geers et al. 
|2006l |2007] l. PAHs are heated by UV-photons from the central 
star. In many cases, the PAH emission is found to be spatially 
extended (e g. van Boekel et al. 2004; Habart et al. 2004 2006) 
Visser et al. |2007l Doucet et al. 2007 & Boersma et al. 2008 ). In 
particular, in Herbig stars, the PAH emission is more extended 
than the continuum (at similar wavelengths) and can be as ex- 
tended as the total size of the disk (Visser et al. I2007I I. In disks 
around less massive stars, given the lower UV flux, the emission 
is weaker and more compact around the central star Three PAH 
emission features are present in the spectral range covered by 
our observations centered respectively at 8.6 /jm, 11.3 fj.m and 
12.7 jim. All the three sources we study here show evidence of 
emission of (at least one of) these features (e.g. van Boekel et al. 
12005 . and figure 1). The PAH emitting region may be much larger 
than the continuum/silicate one (Sect 4.2). As a consequence, at 
the position of the PAH emission there are (at least) two visibil- 
ity components which contribute to the MIDI observed visibil- 
ity. Depending on the relative contribution of PAH emission to 
the total flux at a specific wavelength, the total visibility tends 
to be lower due to the presence of the extended PAH emission. 
Given the contamination of telluric absorption corresponding to 
the 12.7 jum emission band, we do not take this feature into ac- 
count in our analysis. The main findings that emerge from the 
analysis of the MIDI visibility and correlated flux are: 

PAH 8.6 fxm feature 

1. HD 101412: The visibility increases monotonically 
between 8 jim and 9 /im and the disk appears slightly 
more extended shortward of 9 //m than at 10 fim. This 
may be a signature of extended PAH emission 

2. HD 135344 B: No evidence of 8.6 fim emission in the 
total spectrum or in any of the correlated spectra 

3. HD 179218: The visibility is systematically low short- 
ward of 9 fim. However, it is not clear whether this just 
reflects the characteristic sinusoidal modulation of the 
"ring-like" structure or whether it is the PAH emission 
that causes this effect 

PAH 11.3;um feature 

1. HD 101412: A drop in the visibility of HD 101412 
around 11.3 fim is detected along the second baseline 
(59 m, 102°). This is also the observation with the best 
signal-to-noise ratio for HD 101412. A similar trend is 
visible along the third baseline (62 m, 110°), although 
the large error bars make this detection uncertain. HD 
101412 is classified as a group II (self-shadowed) disk. 
The detection of PAH emission more extended than the 
dust emitting region is peculiar for such a disk structure. 
The UV-excited PAH molecules should lie in a dust-free 
surface layer of the disk. This source has a strong [O i] 
6300 A emission - equally unusual for a group II source 
- which led van der Plas et al. ( 120081 ) to consider HD 
101412 to be somewhat transitional in character between 
group I and group II. 
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2. HD 135344 B: A clear drop in visibility in the vicinity 
of 11.3 jjm is detected with MIDI with all baselines. 
This is even more clear in the correlated flux (Fig. |5]l: 
no emission is detected in the four measurements of 
the correlated flux, indicating that the 11.3 fim PAH 
emission comes from spatial scales larger that those 
resolved by MIDI. 

3. HD 179218: No detectable variation of the visibility is 
visible in the MIDI observations of HD 179218 around 
11.3 fim. This disk is classified as a group I source 
(flared). The flaring geometry and the large size found 
for the dust emitting region may suggest that in this case 
dust and PAHs have a similar extent. 



Appendix B: Inclined ring model 

In the text, two simple geometrical models were used to fit the 
MIDI visibility: 1) a uniform disk and 2) a uniform ring. The 
uniform disk model is a first-order approximation of a spherical 
brightness distribution. The intensity (lip, (p), p and (p polar co- 
ordinates) is constant and non zero within a radius 9 12 and null 
outside: 



Hp, 



p < 0/2 
p > 0/2 



given the azimuthal symmetry, I(p, 
form of Eq. IB.ll is: 



(B.l) 



I(p). The Fourier trans- 



V(u, v) = 2 



Ji(7T0r) 



(B.2) 



(ner) 

with Ji Bessel function of the first order. The second model used 
in this paper is the uniform ring. In polar coordinates the bright- 
ness distribution of the uniform ring is: 



f p< 0in/2 

I(p) = \h) 0in/2<p<eou,/2 
p> 0out/2 



(B.3) 



with 0in and 0out inner and outer diameters of the ring. 

Given the linearity property of the Fourier transform, the vis- 
ibility of the uniform ring with inner size (diameter) 0,„ and outer 
size 0out is the diff'erence of the visibility of two uniform disks. 
The visibility of a structure inclined by an angle / and rotated by 
an angle (p (position angle of the major axis) can be obtained by 
rotating the u, v coordinates (e.g. Berger & Segransan l2007l l: 

u' — u ■ cos((f>) + V ■ sin(p (B-4) 
v' — -u ■ sin{(p) + V ■ cos(p (B.5) 

and by applying a compression factor along the minor axis: 

r = ^Ju'^ + v'2 ■ cos{i)^. (B.6) 
The visibility of an inchned uniform ring is: 



Vringir) = /i + 2- 



/2 _ a2 L 

our ^in 



i'r0ou,r) 



(7i0i„r) 



■] (B.7) 



where is the fractional contribution of an unresolved compo- 
nent (e.g. the central star) to the total flux at wavelength A. In 
the extreme case of an infinitely thin ring with 0,„ - 9out, the 
visibility becomes: 



The model adopted for HD 179218 is the sum of two uniform 
rings: 



/(p) = 



^0 


P<0\/2 


k 


0i/2<p< 02/2 





02/2 <p< 0-i/2 


h 


03/2 <p< 04/2 




p > 04/2 



(B.9) 



with 01,02, 03 and 04 the diameters of the two uniform rings re- 
spectively. The visibility of such a model is: 



V(r) = fAVuR.i„{0l,02) + (1 - fA)VuR.ou,(03, 04) 



(B.IO) 



Vn„g(0in) =/)+(!- fA)Mn0i„r). 



(B.8) 



with VuR.„ and VuR„ut the visibility of a uniform ring (Eq. IB. 7b 
for the inner and outer component respectively. 
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